Data import

Clean and mutate features.csv Each row in this file is associated with a track a respondent entered in the mapping tool of the survey (custom programmed in Qualtrics). First, we recode trip frequencies as stated by respondents in form of an interval. We take the mean of the interval for further analysis. Note, respondents stated trip frequencies for the 2018 summer season (current summer) and answered a hypothetical question how the number of trips would have changed would elodea clog some of the waterbodies they boated on.

Calculating annual operating cost per route. Here, we calculate the annual km as the product of the stated trip frequency, round-trip length in km (calculated by ArcGIS), the stated hourly boat operating cost, and the inverse of the assumed average travel speed of 10 miles/h (16km/h). Both pre- and post-invasion costs are calculated.

Missing data, outliers, and recoding on other variables. Associating each route with either private, commercial, or government related operators.

## Warning: NAs introduced by coercion

Preparing features file for leaflet shiny app

Clean and mutate data.csv file

## Warning: NAs introduced by coercion

## Warning: NAs introduced by coercion

Preparing data file for leaflet shiny app

Survey analysis

1. Response

Of the 5095 mailed envelopes, 3.6% were undeliverable for an effective sample size of 4914. Of these the survey received 965 responses for a response rate of 19.6%. The effective sample contained 220 business contacts which resulted in 21 business responses. Some respondents (0.5%) reported to have online difficulties with mapping their boat tracks and 0.3% reported to not have Internet. 21 respondents returned the $1 bill by mailing the bill back to the research team. Two respondents requested to be taken off the mailing list for the reminder mailings.

Of the 965 respondents, 36 percent solely operated in marine waters and consequently were not asked to provide mapping detail of their 2018 boat tracks. 29 percent (n=280) of respondents solely operated in freshwater, while 18 percent (n=175) operated in both fresh and marine waters which explains why the resulting geographic information on boat tracks includes tracks in marine waters. 15 percent reported no to have used their boat in 2018.

Of the 965 respondents, roughly a third (322) gave mapping answers while the remainder largely operated in marine waters (350 of 643 non-map responses) and consequently was not asked about their boating routes. A total of 137 respondents who indicated to have operated in Alaska freshwater did not use the mapping tool. One of the reasons for this non-response to the mapping tool was the tool’s complexity as voiced to the research team by less than 1 percent of responsents.

As expected, respondents with mapping answers took considerable time longer to answer the survey. Response times for these respondents was within the approximated time noted in the consent part of the survey’s introduction section and amounted to 12 minutes if the respondent operated only in freshwater and 14 minutes if operating in fresh and marine waters. Response times for respondents without mapping their routes was 7 minutes for freshwater only and 10 minutes for operators in fresh and marine waters.

Tables: Number of respondents with and without mapping answers

## Warning in kableExtra::kable_styling(.): Please specify format in kable.
## kableExtra can customize either HTML or LaTeX outputs. See https://
## haozhu233.github.io/kableExtra/ for details.
Summary of respondents with mapping response, n=322
User median.min count elodeaCount
Did not operate 9.2 4 0
Fresh and marine 14.0 115 15
Only freshwater 12.4 203 28
Summary of respondents without mapping response, n=643
User median.min count elodeaCount
Did not operate 0.7 156 0
Fresh and marine 9.6 60 6
Only freshwater 6.9 77 11
Only marine 0.9 350 0

1.2 Representativeness of the sample

We conducted an independent t-test to see how representative respondents to the survey were in comparison to the general Alaska population. Ideally we had hoped to compare the sample to the sample frame (population) of registered boaters but there is no numeric variable that could be of use to conduct such a comparison. The only available metric to conduct this comparison is personal income reported in the survey which can be compared to the personal income (PINCP variable) in the PUMS dataset for Alaska associated with the 2017 American Community Survey. As to be expected, the result of the t-test shows dissimilarity in samples (p<0.000, t=47.7, df=2313.3). The mean income is much higher in the sample (91,398) compared to the observed population mean in the ACS (38,624) The Figure below shows the sample distribution of income in red while the ACS income distribution is shown in black.

## 
##  F test to compare two variances
## 
## data:  data6$Income and PUMS$PINCP
## F = 0.31147, num df = 964, denom df = 5159, p-value < 2.2e-16
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
##  0.2830600 0.3438708
## sample estimates:
## ratio of variances 
##           0.311475

## 
##  Welch Two Sample t-test
## 
## data:  data6$Income and PUMS$PINCP
## t = 47.714, df = 2313.3, p-value < 2.2e-16
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  50605.88 54943.87
## sample estimates:
## mean of x mean of y 
##  91398.96  38624.09

1.3 Geographic representativeness

The Qualtrics online survey software collects respondent information such as IP addresses which can be used to determine the location the respondent was in when answering the survey. If the respondent completed the survey using the Qualtrics Offline App on a GPS-enabled device, this data will be an accurate representation of the respondent’s location. For all other respondents, the location is an approximation determined by comparing the participant’s IP address to a location database. Inside the United States, this data is typically accurate to the city level.

Based on this information, we created the following maps. The first map shows that registered boat owners who responded while in Alaska were mainly from the Southcentral Alaska (435), Southeast Alaska (126), Interior (62), Cordova (18), Kenai (26), Kodiak (18), Bristol Bay (7), Y-K Delta (7), Norton Sound (4), and northern Alaska (4).

There was considerable representation from out-of-state responses as well (138) and Hawaii (2). Some of the outside responses may also be by boat owners responding while on vacation or other trips.

Figure: Map showing clustered locations for the 965 respondents

Respondents from lower 48 locations were more likely to report that they did not boat in Alaska in 2018 as shown by the blue markers in the map below. No map response was geographically equally distributed.

Figure: Map of respondents and where they reside color coded by whether repondents gave mapping answers or did not boat (n=848)

2. Elodea and other AIS concern

The survey asked whether respondents took their boats into any elodea infested waterbodies. The most frequented elodea-infested waterbody as reported by 38 respondents was the Chena River (Chena Slough), which consequently is the most concerning elodea infestation subject to boat transmission in the state. Less frequented are the Eyak River (8), Eyak Lake (6), and Alaganik Slough (3). Relatively speaking, these frequencies in the Cordova area are large relative to the population of boaters in this area.

Table: Respondent count for boat owners who took their boat into elodea infested waterbodies in 2018, n=51
Count of respondents visiting waterbodies with known elodea infestations in 2018, n=51
location count
Alaganik Slough 3
Chena Lake 2
Chena River 34
Chena Slough 4
Eyak Lake 6
Eyak River 8
Lake Hood 1
Sand Lake 1
Sports Lake 1
Stormy Lake 3
Totchaket Slough 3

2.1 Elodea transmission risk map

Data on boaters’ mapping responses were used to create the following risk map showing the elodea infestations by size and the reported boat routes. The red circles indicate the relative sive of local elodea infestations as reported in the AKEPIC database areal extent of each infestation. For illustrative purposes, the circle diameter is proportional to log(infested area)*10. This illustration shows how much more extensive the infestations are in the Cordova area compared to infestations in other parts of the state which is consistent with the hypothesis that elodea first established in Cordova.

The map below also shows the extent of routes taken by boaters who are not following a 100% clean drain dry best practice. The extent to not adhering to this best practice is geographically equally distributed across the state. This result suggests that outreach efforts to follow clean drain dry procedures should be broad in scope across Alaska but particularly focused on Fairbanks, the Interior and the Cordova area. Figure: Respondents’ (n=324) boat routes and locations of known elodea occurrence

## 
## Attaching package: 'geojsonio'
## The following object is masked from 'package:jsonlite':
## 
##     validate
## The following object is masked from 'package:base':
## 
##     pretty
## Warning in rgdal::readOGR(input, rgdal::ogrListLayers(input), verbose =
## FALSE, : Dropping null geometries: 156, 555, 684
## Warning in rgdal::readOGR(input, rgdal::ogrListLayers(input), verbose =
## FALSE, : Dropping null geometries: 223, 312

Analysis of the routes taken by respondents who indicated to have boated outside Alaska first before bringing their boat to Alaska freshwater showed that boats were taken to Skilak Lake, Tustumena Lake, the Susitna River out of Talkeetna going North, and the Chena River. There is also one route in marine waters in Port of Valdez and Southeast Alaska.

Figure: Map of boat routes taken by boaters who first boated outside Alaska in 2018, n=4

## Warning in rgdal::readOGR(input, rgdal::ogrListLayers(input), verbose =
## FALSE, : Dropping null geometries: 156, 555, 684

Alaska boat owners’ concerns related to AIS impact on sport fishing are highest followed by biodiversity , subsistence, and recreation values, as shown by the median rank in the next Table. Respondents were less concerned about impacts on businesses, real estate, and commercial fishing as well as boating safety.

Table: Descriptive statistics related to respondents’ rating of concern for AIS affecting specific sectors, n=372
Descriptive statistics related to respondents’ rating of concern for AIS affecting specific sectors, n=372
Statistic Sport fishing Boating safety Recreation value Businesses Real estate values Subsistence Biodiversity Commercial fishing
Min 1 1 1 1 1 1 1 1
Median 5 2 4 3 2 4 4 3
Max 5 5 5 5 5 5 5 5
Pct.Valid 75 64 73 64 68 70 71 66

2.2 Alaska boat owners’ boating and boat purchases occurring outside of Alaska

The survey also contained questions on the purchase of boats from outside of Alaska and whether Alaska boat owners boated outside Alaska before putting their boat into Alaska freshwater in 2018. Only four respondents (less than 1% of respondents reporting, n=311) indicated to have boated outside of Alaska in 2018 before bringing their boat to Alaska. One respondent brought a boat used in Kansas, a state with known zebra mussle (Dreissena polymorpha) infestations. Other outside locations included Washington state, and British Columbia and the Yukon Territories (Table).

Table: Boat owner count for owners who took their boat out of state before boating in Alaska in 2018
Percentage of boat owners boating outside of Alaska in 2018
locOutside freshOutside count percentOfN
No 307 0.9219219
Did not report 59 0.1771772
British Columbia Yes 1 0.0030030
Kansas Yes 1 0.0030030
Washington Yes 1 0.0030030
Yukon Territories Yes 1 0.0030030

Boat owners who take their boats outside of Alaska are more likely to bring AIS into Alaska. We compared the characteristics of boat owners who take their boats outside of Alaska (Table) with the average boat owner responding to the survey (Table). Due to the small sample size of outside boaters, we caution the reader to place too much emphasis on the results. Outside boaters showed lower mean income, are perhaps retired, but operate slightly more expensive boats as shown by the higher than average hourly operating cost. Also, these boat owners had shorter the average annual boating distances.

Table: Boat owner characteristics who boated outside Alaska before boating in Alaska
Summary statistics for operators who took their boats outside Alaska in 2018, n=4
Statistic Income Age Passengers Number of unique routes Km/year Operating cost $/h Operating cost $/year Trips/year Trips/year contingent
Mean 75000 52 2 2 91 25 2590 4 4
Std.Dev 43301 20 1 2 37 11 2141 4 4
CV 1 0 0 1 0 0 1 1 1
Min 37500 29 1 1 53 15 933 2 2
Q1 37500 35 2 1 63 16 1005 2 2
Median 75000 54 2 1 87 22 1946 2 2
Q3 112500 68 2 3 120 34 4175 6 6
Max 112500 70 3 5 138 40 5533 10 10
Pct.Valid 100 100 100 100 100 100 100 100 100

Additional analysis was conducted regarding the number of boats purchased outside Alaska, when and where these boats were purchased and whether they were purchased used or new. Since the survey asked about historical information, this question is most likely biased in terms of whether respondents remembered correctly. This bias may have resulted in underreported purchases that occurred farther in the past and may have contributed to the increase in purchases shown in the Figure below. If we ignore this recall bias, it appears that there is an increased influx of used boats from outside Alaska that comprise significant AIS risk.

Figure: Reported historical purchase of boats from outside Alaska

Boats that are operating in Alaska and were purchased out of state are originating all across the lower 48 with some originating in Canadian provinces and territories. Most concerning are used boats purchased in mussle infested states such as Michigan, Wisconsin, Indiana, Kansas, Texas, and Colorado (Figure). Due to the analysis being based on a survey and not a census of Alaska boat owners, the map below shows only a 2018 snapshot for a selection of boat owners responding to the survey. Likely, used boats are purchased by Alaska’s registered boat owners from all lower 48 states and Canadian provinces and territories.

Map: States where boats were purchased by type of purchase

2.3 Clean Drain Dry

Interesting to note, more than half (61%) of survey respondents who answered the question on whether they follow clean drain dry best practices indicated that they do this all the time. Only 16 percent reported to never follow this practice, and 8% each follow this practice 25 percent, 50 percent, and 75 percent of the time. Figure XA shows that the lowest proportion of 100% compliance with clean drain dry was reported by boaters between the ages of 41 and 50. Figure XB also shows that respondents are normally distributed across ages, with boaters between the ages of 51 and 70 most frequently represented.

Figure: Stated percentage of time respondent followed clean-drain-dry procedure to minimize transmission of aquatic invasive species by age group (n=312).

Figure X is used to look at clean drain dry compliance across age groups and the annual boating distances to see which age group may represent a higher risk category in terms of long-range transmission risk due to not following clean drain dry procedures. Two age groups with longer than average boating distances are the ages between 31 and 40 and ages between 51 and 60. While the former group showed relatively high compliance with clean drain dry, any outreach efforts on clean drain dry would effectively target the age group of 51 to 60 and older.

Figure: Stated percentage of time respondent followed clean-drain-dry procedures to minimize transmission of aquatic invasive species by age and annual distance

## Warning: Removed 37 rows containing missing values (geom_point).

Boater and trip characteristics

Analysis of the recorded boating tracks via the mapping tool of the survey showed that boaters between the ages of 51 and 60 travel the longest annually and show the highest annual trip frequencies. Younger boaters, particularly the ages between 21 and 40, are more likely to have longer boat routes and less frequent trips than older boaters. The weighted average boating distance per trip for the younger age group is 50 to 100% longer than the weighted average boating distance per trip of other age groups, while the annual trip frequency is less than a tenth of the annual trip frequency of the 51-60 age group for example.

Table: Total distance, total trips, total respondents, and weighted average distance per trip by age group in the personal user group
Total distance, trips, respondent count, and weighted average distance per trip by age group
ageBin count sumTotalKm sumTrips wghtAvgDis
21-30 6 4858 117 41.5
31-40 37 20102 576 34.9
41-50 56 15108 678 22.3
51-60 75 35243 1223 28.8
61-70 81 16764 748 22.4
71-80 24 4362 163 26.8
>80 1 925 16 57.8
NA 73 26078 1233 21.2

Not surprisingly the main purpose of all recorded routes in the mapping tool was sport fishing, amounting to half of all 758 routes that were marked by respondents. Non-specified purpose category ranked second with 20 percent of all routes followed by cabin related trips and hunting (9 and 7 percent). Camping is associated with the longest average route distance of 153 km while bird watching and cabin related trips have the lowest average route distances. The purpose that is associated with the highest range of route distances was also camping where routes ranges from 1km to 883 km followed by the other category and sport fishing (0km to 497km).

Table: Statistics related to the main purpose of routes
Summary statistics for the primary purpose of routes, n=324
Primary purpose Route count % of all routes Mean rank Average route in km Shortest route in km Longest route in km
Sport Fishing 384 50 2 42 0 497
Other 149 20 3 40 0 675
Cabin 68 9 1 27 0 139
Hunting 50 7 2 74 3 400
Subsistence Fishing 45 6 2 53 1 478
Camping 25 3 2 153 1 883
Wildlife watching 25 3 2 58 1 257
Hiking 5 1 4 62 18 90
Bird watching 4 1 1 22 0 55
Lodge 2 0 2 71 57 86
Mining 1 0 2 36 36 36

Since the survey was aimed at improved understanding of AIS transmission risk related to Alaska’s boaters, respondents to the mapping tool and their characteristics were of particular interest. Below, we show descriptive statistics to the boaters who plotted their boating routes in the survey’s mapping tool and who boated for personal reasons. Final respondents to the survey is largely comprised of this group of respondents.

As expected, the average respondent’s income (median 87500) is much higher compared to the statewide population average (median 36000). This difference is partly explained by the cost ot operate boats which was estimated for the average boater to equal 3429. Boater age is normally distributed across the sample with a mean of 55. Additionally, the average Alaska registered boat owner operating primarily in fresh water has two unique routes, travels 262km a year, and takes on average 14 trips. The average operating cost per hour is 20/h.

Table: Descriptive statistics related to personal operators who gave mapping answers, n=324
Summary statistics for personal operators with mapping answers, n=324
Statistic Income Age Passengers Number of unique routes Km/year Operating cost $/h Operating cost $/year Trips/year Trips/year contingent
Mean 102590 55 3 2 262 20 3429 14 14
Std.Dev 50369 12 2 2 818 20 7087 25 25
CV 0 0 1 1 3 1 2 2 2
Min 12500 22 0 1 0 0 0 1 1
Q1 62500 45 2 1 14 6 164 2 2
Median 87500 56 2 1 44 13 653 6 6
Q3 137500 65 3 3 190 27 2817 16 16
Max 225000 81 30 11 8036 100 44436 265 265
Pct.Valid 70 79 80 100 100 68 68 100 100

Boat operating cost

Tables: Annual operating cost followed by operating cost per hour by operator type
Estimated annual operating cost by operator type
type count mean median sd max min cv
commercial 7 12848.1 3822.2 17607.9 41060.0 41.0 1.4
government 4 6301.8 6301.8 NaN 6301.8 6301.8 NaN
personal 360 7066.4 968.2 15276.9 113662.1 0.0 2.2
Stated hourly operating cost by operator type
type count mean median sd max min cv
commercial 8 21.2 22.2 18.9 40 1.5 0.9
government 15 29.4 5.0 33.6 70 5.0 1.1
personal 741 21.0 15.0 18.5 100 0.0 0.9